Structure, elastic moduli and thermodynamics of sodium and potassium at ultra-high 

pressures 
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The equations of state at room temperature as well as 
the energies of crystal structures up to pressures exceeding 
100 GPa are calculated for Na and K . It is shown that the 
allowance for generalized gradient corrections (GGA) in the 
density functional method provides a precision description of 
the equation of state for Na, which can be used for the calibra- 
tion of pressure scale. It is established that the close-packed 
structures and BCC structure are not energetically advanta- 
geous at high enough compressions. Sharply non-monotonous 
pressure dependences of elastic moduli for Na and K are pre- 
dicted and melting temperatures at high pressures are esti- 
mated from various melting criteria. The phase diagram of K 
is calculated and found to be in good agreement with experi- 
ment. 
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The theoretical and experimental studies of the matter 
properties at ultra-high pressures arouse a great interest 
in the connection with the possibility to obtain phases 
with uncommon properties as well as geophysical and as- 
trophysical applications. As an example, the problem of 
metallic hydrogen can be mentioned!!]. In the high pres- 
sure studies the alkali metals can be conveniently used as 
model objects. This is due, first, to their high compress- 
ibility and, second, to the variety of physical phenomena 
occurring in their compression and numerous structural 
and electron phase transitions (see, e.g.B~@). For heavy 
alkali metals it is th&famous s—d isostructural FCC-FCC 
transition (see, e.g. ,E£l and references therein) as well as 
the transitions to uncommon distorted phases at higher 
pressurestil. Recently it was supposed, basing on the elec- 
tron structure calculations, that lithium can transform at 
high enough-aressures into "exotic" phases similar to that 
of hydroger£3. Thus, further theoretical investigations of 
structural properties of alkali metals at ultra-high pres- 
sures seem to be interesting and important. 

Despite a lot of considerations, this is still an open 
problem. The most of early attempts used computa- 
tional approaches which were not accurate enough from 
the contemporary point of view. It is well known (see, 
e.g.jEJ) that the highly accurate quantitative description 
of the electronic and, especially, lattice properties of met- 
als needs the consideration of the real form of potential 
in the crystal and going beyond the frame of local ap- 
proximation in the density functional, in particular, the 



allowance for generalized gradient corrections (GGA)0. 
In the present work a consistent theoretical study of the 
relative stability of crystal structures of Na and K under 
pressure as well as a variety of related lattice properties, 
is performed basing on these first-principle calculations. 
The most interesting result obtained is that, contrary 
to the traditional concepts (see, e.g.Ja) neither structure, 
which is characteristic of metals under normal conditions 
(BCC, FCC and HCP), is stable at high enough pressures 
even in Na where there are no electron transitions. 

The ah initio calculations of electronic structure, ther- 
modynamical potential, equilibrium lattice parameters 
and elastic moduli at temperature_T = were carried 
out using the FP-LMTO methodE^I with allowance for 
the GGA in the form proposed inta. A-,careful optimiza- 
tion of the parameters of this methodEj made it possible 
to carry out the calculations of the total energy with an 
accuracy within the limits of 0.1 mRy/atom. Parameter 
c/a for the HCP lattice was determined by the minimiza- 
tion of the total energy for a fixed specific volume, and 
the elastic moduli — by numerical differentiation of the 
total energy with respect to tetragonal and trigonal de- 
formations (see, e.g.Ja). Up to now, there are only few 
works devoted to the first-principle calculations ofi-alas- 
tic moduli of metals under high pressures (see, e.g.,E3 for 
Mo and W). The expressions for the free energy F (con- 
nected with Gibbs thermodynamical potential G by the 
relationship G = F + PV, P = -dF/dV, where P is 
the pressure, V is the volume) and elastic moduli Cap^s 
at finite temperatures can be presented in the following 
form: 



F(V,T) = E ee (V) + F ph (V,T) 



(1) 
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is the free energy 



of the phonon subsystem in the harmonic approxima- 
tion, E e e(V) is the total energy of the electrop, subsys- 
tem obtained from the FP-LMTO calculationslij, P ph = 
—dFph/dV is the phonon pressure, Bq is the bulk mod- 
ulus of the electron subsystem at T = 0, £, q, u> are the 
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number of phonon branch, wave vector and phonon fre- 
quency, respectively. In the expression for the elastic 
moduli C a pys(V, T) the first term corresponds to the elec- 
tron contribution at T = 0, second one - to the quasihar- 
monic contribution due to the effects of thermal expan- 
sion, and third one - to the phonon contribution obtained 
from the differentiation of F p h with respect to the corre- 
sponding deformation parameters. For the calculation of 
phonon contributions to thermodynamical functions the 
pseudopotential model described inliZI, which describes 
with a high accuracy a wide range of lattice properties of 
alkali metals, was used. 

Figs. 1,2 show the results of calculations of Gibbs po- 
tentials at T = O for the BCC, FCC and HCP phases of 
sodium and potassium, respectively. It should be noted 
that in the case of Na the phonon contribution to AG 
(contribution of zero-point vibrations, AG zp ) are G^ zp c — 
G b z c p c = 1.31 x 1CT 5 Ry/atom, G^p - G b z c p c = 1.35 x 1CT 5 
Ry/atom. This is well comparable with the electron con- 
tributions to AG under the normal conditions. How- 
ever, already at P > 1 GPa for Na and practically at 
all the pressures for K the contribution of zero-point 
oscillations to AG can be neglected. Generally speak- 
ing, energy differences of order of 10 -5 Ry/atom is too 
small to be accurately derived in our first-principle cal- 
culations; nevertheless, we have obtained correct phase 
diagram even for sodium at low pressures. In accor- 
dance with the results of calculations, at P = the 
BCC phase and HCP phase have the lowest energy for 
K and Na, respectively (actually, under these conditions 
Na has not HCP but 9R structure whose energy, how- 
ever, is very close to HCP0). It is important to empha- 
size that this difference of Na from K is purely quanti- 
tative: according to the results, shown in the insert to 
Fig. 2, K would have to transit to the hexagonal close- 
packed phase at the negative pressure of the order of 
several kilobars. As the calculations show potassium, un- 
like sodium, transits from the BCC to FCC structure at 
P w 11.6 GPa, which is in excellent agreement with the 
experimental datatffl. In this case the relative change in 
the volume A = (Vb cc — Vf cc )/Vb cc w 0.0067 takes place 
at Vo/V = 2.14. Here and below V) is the experimen- 
tal value of the specific volume at the atmospheric pres- 
sure and temperature 10K, equal to 484.12 a.u.ES. Our 
calculations make it possible to suppose that the differ- 
ence between Na and K is associated with the electronic 
topological transition occurring in the BCC potassium at 
Vo/V k 2 and destabilizing the BCC structure. The sim- 
ilar situation takes place in LiEl while in the BCC sodium, 
within the whole range of pressures, the Van-Hove singu- 
larity goes away from the Fermi level under the compres- 
sion. Generally, sodium seems to be a unique metal in the 
Periodical Table: in the whole region of the existence of 
BCC structure it has no singularities of electronic struc- 
ture near the Fermi level, and the Fermi surface remains 
approximately spherical. 

The calculation results of equations of state for sodium 



and potassium are shown in Fig. 3 along with the exper- 
imental data available. It should be pointed out that 
the experimental data agree with the theoretical results 
within their accuracy limits (« 10%). This creates the 
prerequisites for the development of pressure scale based 
on sodium as a reference substance. Note also that at 
room temperature the role of the phonon contribution 
to pressure falls under the compression, and this contri- 
bution in itself is small (of the order of 0.3 GPa at full 
pressure 20-30 GPa). 

Figs. 4,5 display the calculation results of the depen- 
dence of elastic moduli Cij(V) on the compression for 
sodium and potassium, respectively. A drastically non- 
monotonous behavior of shear moduli associated with the 
tetragonal (C — (Gn — C\2)/2) and trigonal (G44) de- 
formations in both BCC and FCC structures is notice- 
able. It should be pointed out that at least at compres- 
sions Vo/V < 2 the calculation results of Cik(V) in the 
pseudopotential model and in the first-principle approach 
are close. In this region the equations of state coincide 
in these two approaches with an accuracy up to several 
percents. This confirms a sufficiently high reliability of 
our use of pseudopotential model for the calculation of 
phonon contributions to thermodynamical values. Nev- 
ertheless, the phonon contributions to the shear moduli 
do not exceed 10% within the whole pressure range stud- 
ied. It should be noted that softening of modulus C is 
a typical pre-transition phenomenon connected with the 
structural-transitions between the BCC and close-packed 
structurealj. However, the softening of modulus G44 in 
the FCC structure of K at high pressures (Fig. 5) is rather 
surprising. It appears to be similar to the softening of this 
modulus, taking place in the_FCC structure of Cs near 
the electron s — » d transitions and is due to the crawling 
of the Fermi level over the peak of d-state density. 

Fig. 6 shows an experimental phase diagram of potas- 
sium and the phase diagram built on the basis of our cal- 
culations. The dependence of the melting temperature 
on pressure, T m (P) in the BCC and FCC phases was 
obtained using the phonon spectra and different melting 
criteria. First of all, we use the Lindeman criterion 

x 2 (T m )/d 2 = const, (3) 

where x 2 (T) — 2M7 - cot ^ ~5t^ * s mean square 

of atom displacement, is the polarization vector, M 
is the atom mass, d is the distance between the nearest 
neighbors. Although the Lindeman criterion is empirical, 
it may be expected that its use for finding the melting 
temperature at high jaressures would be as successful as 
at low temperaturescJ. Nevertheless one can see from 
Fig. 6 that it is not too accuratp-jn a broad pressure re- 
gion. Varshni melting criterion^ which is based on the 
temperature softening of the shear moduli, namely 

G 44 (T m ) /G 44 (0) = 0.65, (4) 

appears to be much more accurate. Here we use the 
method of the calculation of the temperature dependence 
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of elastic moduli from the phonon spectra described inH. 
Note also that we describe with high accuracy the BCC- 
FCC phase boundary. We also present the results ob- 
tained in generalized Debye modeled when all the thermo- 
dynamical quantities are calculated in the Debye model 
but with the Debye temperature found from ab initio 
elastic moduli. One can see that this description is also 
rather accurate for potassium. 

Fig.l shows that the BCC phase of sodium becomes en- 
ergetically unfavorable as compared with the HCP at a 
pressure about 80 GPa. At P > 100 GPa, however, this 
phase demonstrates anomalies in the equilibrium value 
of parameter c/a (Fig. 7). A sharp decrease in the ra- 
tio c/a to 1.2-1.3 at Vq/V > 4.35, which is necessary to 
maintain the HCP lattice in equilibrium, is doubtful ac- 
tually, and seems to be indicative of transition to some 
non closed packed phase with a large number of atoms 
per cell. These .aliases are observed in K, Rb and Cs at 
high pressurestmEil. It is usual to associate their appear- 
ance in heavy alkali metals with the s — > d transition. 
Thus, according to our results, all the three "typically 
metallic" structures - BCC, FCC, HCP do not corre- 
spond to the lowest energy in Na where no electronic 
transitions are observed within the pressure range con- 
sidered. In order to understand qualitatively the cause of 
appearance of " nonstandard" metallic phases, let us use 
the above mentioned pseudopotential model for the esti- 
mations. In this model the radius of "hard" ion core is 
described by the pseudopotential parameter rji3. As the 
estimates show the compression Vq/V ~ 4 correspond- 
ing to the instability of the close- packed phase coincides 
with the condition of overlapping ion cores 2ro ~ d for 
Na. Hence, the concept of well determined ion cores, be- 
ing at the base of standard metallic bond description, be- 
comes inapplicable at ultra- high pressures. As a result, 
as we have seen, the substance transforms into exotic 
non close-packed phases. These results are in qualitative 
agreement with the resultala for lithium. 

In conclusion, note that it would be interesting to 
study the structure of sodium at ultra-high pressures, 
which, as follows from the results obtained, may prove 
to be surprising. Another result of this work, permitting 
a direct experimental check is non-monotonous behavior 
of Na and K shear moduli at pressure. At last, pre- 
cision theoretical description of the equation of state of 
sodium would make possible to use it for the development 
of an accurate pressure scale up to 100 GPa. Although 
the contemporary first-principle calculations can provide 
high enough accuracy also for another substances (see, 
e.g., recent calculations^!! for Si) a very high compress- 
ibility of sodium and the absence of phase transitions in 
a broad range of pressures make it probably the most 
suitable for these purposes. 

The authors are grateful to D. Yu. Savrasov and S. Yu. 
Savrasov for the permission to use the author's version 
of the code realizing the methodt3 in their work as well 
as to D. Yu. Savrasov and E. G. Maksimov for useful 
discussions of the details of this method. 
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Figure captions 

Fig.l. Pressure dependence of the differences of Gibbs 
potentials between BCC and FCC as well as HCP and 
FCC structures for Na. 

Fig. 2. Pressure dependence of the differences of Gibbs 
potentials between BCC and FCC as well as HCP and 
FCC structures for K. 

Fig. 3. Equations of states for sodium and potas- 
sium at T = 295|Kn. Solid line corresponds to FP- 
LMTO calculationsta, dashed liee - to the calculations 
by the pseudopotential methodliU, Empty (solid) trian- 
gles ,|-cii|cles and asterisks (squares) are the experimental 
dataE5c4El for Na and K, respectively. 

Fig. 4. The dependence of elastic moduli C and C44 
on the compression U for BCC Na; empty circles and 
squares show, respectively, the data fromtl; the solid ones 
- the values obtained in the present work. 

Fig. 5. The dependence of elastic moduli C and C44 on 
the compression U for BCC and FCC phases of K. Solid 
(empty) circles and squares show C' and C44 values for 
BCC (FCC) phases, respectively^The dashed line — the 
C values for the bcc phase fromtl 

Fig. 6. Phase diagram of potassium. The solid line - 
experimental dataQ, dashed line- the calculations using 
Varshni criterion (0) dashed-dot line - the calculations 
using Lindeman criterion (^), dotted line - generalized 
Debye model (see the text). Solid circles - BCC-FCC 
phase boundary from our calculations. 

Fig. 7. Dependence of the total energy of HCP struc- 
ture for Na on the ratio c/a for various compressions: the 
solid line — U = 0.75; dashed line — U — 0.76; dashed- 
dot line — U = 0.765. The insert shows the equilibirum 
values of the parameter c/a for HCP structure depending 
on U is shown. Solid (empty) circles denote the values 
taken in the global (local) minimum of the total energy, 
correspondingly. 
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